clear
n=1;
for j=0.00:0.2:8
    results(n,:)=[j two_regimes(j)];
    n=n+1;
end

xx=results(:,1);
yy1=results(:,2);
yy2=results(:,3);
yy3=results(:,4);
yy4=results(:,5);


figure
FSsize=14;
plot(xx+3,yy3,xx+3,yy4,'Linewidth',3)
box('off')
xlabel('Coefficient of variation within volatility regime 2, $\theta_2$','Interpreter','latex','FontSize',FSsize)
ylabel('$X_1^\ast, X_2^\ast $','Interpreter','latex','FontSize',FSsize,'Rotation',0)

annotation('textarrow',[0.5 0.5],[0.6 0.6], 'String','$X_1^\ast$','Interpreter','latex','FontSize',FSsize)
annotation('textarrow',[0.5 0.5],[0.6 0.6], 'String','$X_2^\ast$','Interpreter','latex','FontSize',FSsize)

